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PSEUDO-ABELIAN INTEGRALS ON SLOW-FAST 
DARBOUX SYSTEMS 

MARGIN BOBIENSKI, PAVAO MARDESIC, AND DMITRY NOVIKOV 



Abstract. We study pseudo-Abelian integrals associated with 
polynomial deformations of slow-fast Darboux integrable systems. 
Under some assumptions we prove local boundedness of the num- 
ber of their zeros. 



1. Introduction and main result 



Pseudo Abelian integrals appear as the principal (linear) part of 

,_, the displacement function in polynomial deformations of Darboux in- 

C/5 tegrable cases. This paper is a part of the program of proving uni- 

Q form finiteness of the number of zeros of pseudo Abelian integrals, see 

4^ IH IBM[ IBMNj . After studying the generic cases [Nl IBMj , nongeneric 

C^ cases must be studied, too. Here we study zeros of pseudo Abelian 

Q integrals associated to deformations of slow-fast Darboux integrable 

""^ systems. 

i-H More precisely consider Darboux integrable system uoq given by 

S -o = Mf^, (1) 

O ' 

CA where 

• 

o 
o 



M = l[Pi, Ho = l[P^% ai>0, Piem[x,y]. (2) 



^ We consider the family of forms u^ given by 

•i-H 

H ^^ = p^M^ + eMdPo, PoeR[x,y]. (3) 

Note that for ^ = this form defines a foliation singular along the 
curve Po = 0. The form ^ is Darboux integrable with first integral 

H, = HoP^. (4) 

The system (|3]) is a slow- fast and Pq = is the slow manifold. The 
fast dynamics is given by the Darboux integrable system (fTl). Figure 
|1.1 1 presents typical example that we treat. 



This research was supported by ANR ANAR Project, Polish MNiSzW Grant No 
N N201 397937, and by Israel Science Foundation (grant No. 1501/08). 
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Figure 1.1. Phase portrait of the integrable system. 

Assume that the system ^ has a family ^e{t) C H~^{t) of cycles. 
Consider the polynomial perturbation of the system ^ . 

U^^S = ^E + eiT], £i > 0. (5) 

The linearization in perturbation parameter Si of the Poincare first 
return map is given by the pseudo-Abelian integral 

m=f ^- (6) 

In this paper we study pseudo-Abelian integrals taken along the cy- 
cles ^e{t) of the simplest integrable system bifurcating from a slow- fast 
system: 

Assume that a compact region D is bounded by Pq = and some sep- 
aratrices P^ = 0, i = 1, . . . , A;. Assume that the functions P^, i = 0, . . . A: 
are smooth and intersect transversally in D and that the foliation ([l]) 
has no singularities on IntD. 

Assume moreover that Pq = is transverse to the foliation cjq = in 
all points of P^ n {Pq = 0} except for one point po, where the contact 
is quadratic. Then for ^ 7^ a singular point p^ bifurcates from pq. It 
corresponds to a real value t^ = H^{p^). The bifurcating singular point 
p^ is a center entering the domain P for 6: > 0. Let 7£(t) be the family 
of cycles in the basin of the center bifurcating from p^. For each 5, the 
cycles are defined on an interval [0,^^]. 

Remark 1. Note that if D is compact, then necessarily {Pq = 0} must 
have at least one contact point with the foliation cjq = 0. Below in the 
paper we call turning point the point of tangency between {Pq = 0} 
and the leaves of cjq = 0. 

For ^ < instead of a center, a saddle point bifurcates from p^ just 
outside of D. 

Example 1. Consider the foliation 

uj, = PodPi + ^PidPo, Po = ?/ - x^ Pi = 1 - ?/, (7) 
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with the first integral H^ = {1 — y)[y — x^Y . It has a critical point 
Pe = (1TJ5 0), which is a saddle for 6: < 0, coincides with the tangency 
point po = (0, 0) for ^ = 0, and is a center for 6: > 0. For each s > 
the region D = {x^ < y < 1} C M^ is filled by a nest of cycles %{t) 
vanishing at p^. These cycles are parameterized by values of t varying 
from (corresponds to the polycycle dD) to 



te = H,{p,) = {l + e)-Ul + 



1 



-^ e 



as 6: ^ 0. 
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Figure 1.2. Phase portrait of i7^ = (1 - y){y - x^Y 



Theorem 1. Let I^{t) be the family of pseudo-Abelian integrals as de- 
fined above with the above genericity assumptions. Then there exists a 
bound for the number of isolated zeros of the pseudo-Abelian integrals 
I{t^e), for £ > O and t G (0,^^). The bound is locally uniform with 
respect to all parameters: the form varying in some finite dimensional 
analytic family^ the coefficients of the polynomials^ exponents and in 
particular with respect to e. 

Remark 2. The estimates on the number of zeros of pseudo-Abelian 
integrals do not imply in general estimates on cyclicity of correspond- 
ing polycycles, due to presence of the so-called "alien" cycles, see e.g. 

In order to prove the theorem we want to apply variation and ar- 
gument principle as in ^BMi iBi iBMN j . The difficulty lies in the fact 
that the family of cycles of integration ^^{t)^ e > Q does not have a 
smooth limit as 6: ^ 0. The family tends to a family of broken cycles 
formed by parts of the slow manifold and parts of leaves of the fast 
foliation. After taking an 6:-scaled variation, the cycle of integration 
^e{t) is replaced by a figure eight cycle Se{t)^ see Proposition [l| and 
the integral I^{t) by the integral J^{t) taken along d^{t). For t not too 
close to the center value t^, the cycle 5^{t) can be moved along leaves 
of the foliation to be at some positive 6:-independent distance from the 
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slow manifold. This property, which does not hold for the initial cycle 
7£(t), implies good analytic properties of J . 

However, as t ^ t^, and ^ ^ 0, both cycles 7£(t) and 5^{t) approach 
the turning point. To overcome this difficulty in a 6:-uniform way, we 
perform a blow-up in the family uj^. To each chart of the blow up 
corresponds a time scale. In blown-up coordinates the cycles 7£(t) and 
5^{t) have smooth limits in respective charts. This proves analyticity 
of / or J in the timescale of the convenient chart. 

1.1. Acknowledgment. The authors express their gratitude to Insti- 
tute of Mathematics of Warsaw University and Weizmann Institute of 
Science for hospitality. 

2. 6:- VARIATION OF THE 7^(t) CYCLE 

The key point of our approach is to study analytic properties of the 
integral along the family of figure eight loops. 

On the smooth leaf Hq = t of the unperturbed equation there exists 
a figure eight cycle lying close to the real segment {Hq = t} H D and 
winding around the points of intersection of this leaf with the curve 



Po = 0, see Figure |2.1[ Since this cycle lies on a finite distance from 
this curve, i.e. in a domain where the leaves of the foliation H^ = t 
depend analytically on 6:,t, one can move it to a complex leaf of the 
foliation H^ = t^ uniquely up to a small homotopy. The figure eight 
family of cycles 6^{t) is defined as the family of these lifts. 



Figure 2.1. The "figure 8" loop 

Let ^e{t) C H~^{t) as above and I^{t) be the integral ^ along ^e{t)' 

Proposition 1. The e-variation of the cycle ^e{t) is equal to the figure 
eight cycle S^{t): 

Var,(7.(t)) = %{t e-^) - 7.(t e— ^) = 5,{t). (8) 

Proof. As in |BM[ IBMNj . we transport the cycle %{t) for t h^ t e""'^ 
and subtract the transport along t \-^ t e~^^^. Using the local expression 
of H^ near the slow manifold, we see that the parts of %{t e^^^) and 
j^{t e~^^^) along the slow manifold cancel. The figure eight cycle 5^{t) 
remains. D 
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Corollary 1. The s -variation of the pseudo-Ahelian integral I^{t) is 
an integral of the form -p^ along the figure eight cycle S^{t). 

{Varjm : = /,(t e^ " ^t e'^ = / ^ = Mt). (9) 

3. Blowing up the turning point 

We need to prove analytic properties of the integrals / and J in a 
neighborhood of the slow manifold with respect to a convenient time 
scale. This will be achieved by a convenient blow-up of the family in 
a neighborhood of the turning point. This section is dedicated to this 
geometric construction. 

Note that by Morse lemma we can put our family ([3]) to a normal 
form of Example [T] in a neighborhood of the turning point: 

Lemma 1. By an analytic e-independent change of coordinates defined 
in a neighborhood of the turning point po we can assume that cjq = dy 
and Pq = y — x'^. 

We want to study the analytical properties of the foliation in a neigh- 
borhood of the turning point (the vertex of parabola) . The difficulty 
is the approaching of the center to the vertex of the parabola, linearly 
with £. So we make a blow-up of the of the turning point of the family 
in the product space {x^y^e) of phase and parameter spaces. The fam- 
ily blow-ups were studied by Dumortier and Roussarie e.g. in [DRJ. 
This is needed since we want to prove analyticity with respect to both 
phase and parameter values. We want our blow-up to preserve the 
parabola y = x^ and to separate the newborn center from the vertex. 
This requirements lead to the weighted blow-up with weights (1:2:2). 

Recall the construction of this weighted blow-up. We define the 
weighted projective space CPf.2.2 as the factor space of C^ by the C* 
action [x^y^e] h^ {tx^t^y^t^e). The blow-up of C^ at the origin is 
defined as the incidence three dimensional manifold M = {{p^q)\q G 
P} C CPl,,, X C\ 

Note that the weighted action of C* on C^ \ {0} is not free. The 
stabilizer of points (0, ?/, s) is a Z/2 = {±1} subgroup. As a result, the 
quotient is not smooth on the line {x = 0}. In chart described below 
we will work with the double covering which is smooth. 

The blow-down tt : M ^ C^ is just the restriction to M of the 
projection CPf^2:2 x C^ ^ C^ 

For future applications we will need explicit formulae for the blow- 
up in the standard affine charts of M. The projective space CPf.2.2 
is covered by three affine charts: Ui = {x j^ 0} with coordinates 
(Yi.Ei), U2 = {y ^ 0} with coordinates (^2,^2) and Us = {s ^ 0} 
with coordinates (Xs^Y^). The transition formulae follow from the 
requirement that the points (l,Yi,£'i), (^2,1, £"2) and (X3, 1^3,1) lie 
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on the same orbit of the action: 



f ,^' (10) 

Remark 3. Note that transition functions are singular since we cover 
smooth double covering instead of the exceptional divisor itself. It is 
easy to observe that replacing coordinates X2 and X3 in charts U2 and 
Us respectively, we obtain usual, smooth projective transition func- 
tions. 

These afhne charts define affine charts on M, with coordinates (Yi, £"1, ti), 
(X2,£'2,t2) and (X3, 1^3,^3). The projection n is written as 

x = ti y = tlYi e = tlEi (11) 

x = t2X2 y = tl s = tlE2 (12) 

x = hXs y = tlYs e = tl (13) 

We apply this blow-up n to the one-dimensional foliation J^ on C^ 
given by the intersection of rf^ = and u^ = 0. This foliation has a 
complicated singularity at the origin. Denote by i{~^T the lifting of 
the foliation T to the complement of the exceptional divisor 7r~^(0). 
This foliation is regular outside of the preimage of the parabola V = 

Proposition 2. The foliation i{~^T can he extended analytically to 
the exceptional divisor 7t~^{0). The resulting foliation 7r*J^ is regular 
outside of the strict transform of the parabola V , the family of centers 
{X3 = 0,13 = Y+ts} ^^^ ^ point ^5 = (1 : : 0) on the exceptional 
divisor. 

Remark 4. The additional singular point has clear geometric inter- 
pretation which is characteristic to the family blow-up of a slow-fast 
system. It is a "trace" of the fast direction on the blow-up of the slow 
manifold. Consider the following toy example of foliation in M^ defined 
by 1-form uj^ = xd(l — y + ax) + £{1 — y + ax)dx. Making the family 
blow up 

tt: y = xY^ £ = x E 

we obtain a foliation given by 2-form 

^7r%uj,Ad£) = {-a-E-aEx+Y+ExY)dEAdx+xdEAdY+EdxAdY 

that vanishes at the point x = 0^ E = O^Y = a corresponding to the 
direction of the fast system. 
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Figure 3.1. The foliation t{~^T 



Proof. We check it in each chart separately. Note that a codimension 
2 foliation J^ in 3-dimensional space is uniquely defined by 2-form a by 
the condition 

J^ = {v : iya = 0}. 

Alternatively, locally a = iy Vol, where Vol is a suitable non- vanishing 
volume form. Singular points of foliation correspond to zeros of 2-form 
a. The foliation ([7| in 3-dimensional space C^ is given by a^ = uj^Ads. 
The pull-back foliation tt* J^ (strict transform of J^) is defined by the 
pull-back 7r*cr£ divided by a suitable power of the function defining the 
exceptional divisor. In charts [/j, j = 1, 2, 3 we have 7r*cr = tj5j, where 

ai =2{Yi - l){Yi -Ei + EitlYi) dE^ A dti+ 

+ {Yi-l-Ei + EitlYi) (tidEi A dYi + 2Eidti A dY^ , 

G2 =2(1 - Xl){l -E2 + E^tl) dE2 A dt2+ (14) 

+ 2(1 - tl)X2E2 (t2dE2 A dX2 + 2E2dt2 A dX^ , 

53 =4X3(1 - t^y3)dt3 A dX3 - 2(1 + Xl -Y^- tlY^)dh A dy3 

The zero locus of the form a in a neighborhood of the exceptional 
divisor consists of germs of two curves and the singular point p^ = [1 : 
: 0] generated by weighted action. These curves are {X2 = ±1, £"2 = 
0) (strict transform of the parabola P = 0) and (X3 = 0, I3 

(family of centers) . 

In chart U'^ we observe that the lifted foliation is given by two first 
integrals: ^3 and 

s = {l-tlY,Y''hY,-Xl) 



i+tl 
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which can be analytically continued to ta = 0. This foliation has no 
singularities near the exceptional divisor except for the line of centers 
X3 = 0,13 = TX72. Note that the strict transform of the parabola is 
outside of this chart. 

D 

4. Proof of the Theorem 

In this section we first take benefit from the blowing-up in the family 
performed in the previous section to prove analyticity of the integrals 
/ and J in convenient time scale. Let us be given a compact family of 
cycles on the exceptional divisor of the blown-up foliation at a finite 
distance from the singularities. We can extend it to a continuous fam- 
ily in its full neighborhood in the total blown-up space. By analyticity 
of blown-up foliation, see Proposition [2} the integrals along these cy- 
cles will depend analytically on the cycle in the extended family. In 
particular, this means that: 

Lemma 2. For any Rg < 5c(0) = e~^ 

(1) the integral I is an analytic function of s and t^ in some neigh- 
borhood of the arc {0}x{Rse'^'^\\(p\ < tt} in complex (^3, s)-plane; 

(2) the integral I is an analytic function of s and ts in some neigh- 
borhood of the segment {0} x [i?s, Sc{0)] in complex (^3, s)-plane; 

Proof of Lemma \^ Note that the restriction of the foliation tt* J^ to the 
exceptional divisor has first integral s = e~^^ {Ys — X|). It is easy to 
construct a compact family of cycles 7 lying on the exceptional divisor 
and corresponding to the values of s mentioned in the first claim of the 
Lemma, so above argument proves the first claim. 

For the second claim consider the compact family of real cycles van- 
ishing at the center. Analyticity of the integral along these cycles is 
standard. D 

Consider the transversal X2 = to the line ^2 = ^2 = in chart 
U2. To each point on this transversal corresponds a figure eight cycle 
as in Proposition [l] passing through this point and lying on a leaf of 
the foliation 7r*J^ (note that this line is a leaf of 7r*J^), so the integral 
J defines a function on this transversal in a neighborhood of 0. 

Lemma 3. The function J is an analytic function 0/^2, E2 for |£'2|, |t2| ^ 
p for some sufficiently small p > 0. 



Proof On the line t2 = E2 = lying in U2 the blown-up foliation (14) 



has two singular points X2 = ±1, see Figure [3Tj The figure eight cycle 
lying on this line is on finite distance from these two singularities. By 
integrability of 7r*J^, it can be extended to all sufficiently close leaves, 
forming a continuous family of figure eight cycles. 

n 
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Proof of Theorem [I[ The first integral maps tlie open nest of vanisliing 
cycles to the interval (0,^^). We split the interval into two parts, one 
being the image /q of the interval [i?5,5c), and the remaining part. 
On the first part the number of zeros of / is uniformly bounded by 
Gabrielov's theorem, due to Lemma [2] (2). 

To estimate the number of zeros on the remaining part we apply 
argument principle to the contour F consisting of two arcs A^ and a^ 



and two straight segments L± joining their ends, see Figure 4.F The 
arc a^ of angle 2s has infinitesimally small radius and the arc A^ is the 
image of the arc described in Lemma [2J(1). 

The ^-uniform bound for the increment of argument along A^ is a 
direct consequence of the Lemma [2| (2). 

The increment of argument along segments L± can be estimated by 
number of zeros of the imaginary part of the function I on L±. We 
calculate using fact that function / is real on the real segment: 



-AArgLJ < #{t : Iml = 0} + 1 = #{t : / - / = 0} + 1 = 

TT 

#{t : I{e^'' t) - I{e-^'' t) = 0} + 1 = #{t : J{t) = 0} + L 



Thus, one translates problem into estimating the number of zeros of 
the integral J along the family of figure eight cycles corresponding to 
segments L±. 

We split the segment L± into two parts. One, closer to tc, is in 
the image of the transversal considered in Lemma (S) The functions J 
on this part can be considered as an analytic function of ^2,^2- The 
remaining part of L± corresponds to the part of the family of figure 
eight cycles lying not closer than some fixed positive distance from 
the vertex of the parabola. These cycles can be deformed along leaves 
to be not closer than some positive constant from the parabola, see 
Figure [2TJ Therefore the function J is an analytic in t, e on this part. 

By Gabrielov theorem, the number of zeros of J on L± is uniformly 
bounded in e. 

To obtain an upper bound for increment of argument of / along the 
small arc a^ we use an argument similar to those from [BM*. ' BMNl iBl IN] . 
One can easily prove that |/| < C\t\^/^ in sectors. On the other hand, 
for any fixed e the function / has a leading term at t = of the form 
t^ log-^ t, see the above papers. Together it proves the existence of the 
uniform in e upper bound for the increment of argument. 

All the above constructions depend analytically on parameters like 
coefficients of the polynomials P^, exponents a^ and coefficients of the 
form rj. 

n 
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Figure 4.1. Contour T 

5. Concluding remarks and open problems 

Remark 5. Consider the model Example [T] An alternative way to 
prove Theorem [T] for this system would be to reduce it to the situation 
considered in our previous paper [BMN]. Indeed, change of variables 
X = y/sX^ y = sY the system Q has the first integral s = {Y — X'^){1 — 
sVy/^. After additional blow-up of the point at infinity, we get an 
unfolding of a system with two saddle-nodes. Unfolding of systems with 
polycycle with two saddle-nodes was considered in |BMNj . However, 
this approach does not generalizes to more general systems considered 
in Theorem [H 

In the spirit of our program of proving of uniform finiteness of the 
number of zeros of pseudo-Abelian integrals, more general slow- fast 
Darboux systems should be studied. Already in the generic, stable 
under small perturbations of coefficients cases one has to face several 
problems, in particular: 

(1) The unperturbed Darboux system cjq = can have extra sin- 
gular points inside D^ not lying on the zero level of Hq] 

(2) Unperturbed Darboux system can have a nest of cycles, and 
the slow manifold cuts the nest and becomes a part of polycycle 
bounding the new nest; 

(3) The curve {Pq = 0} can have additional tangency points with 
the leaves of cjq = 0, which generate saddles type singularities; 

(4) The nest of cycles accumulating to the polycycle and encircling 
more than one newborn singular point, so called "big cycles"; 



These possible scenarios are illustrated on Figure [5T] below, keeping 
numeration. 
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Figure 5.1. Phase portrait of the system with two factors. 
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